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Nonlinear coherent modes are the collective states of trapped Bose atoms, corresponding to difler- 
ent energy levels. These modes can be created starting from the ground state condensate that can 
be excited by means of a resonant alternating field. A thorough theory for the resonant excitation 
of the coherent modes is presented. The necessary and sufficient conditions for the feasibility of 
this process are found. Temporal behaviour of fractional populations and of relative phases exhibits 
dynamic critical phenomena on a critical line of the parametric manifold. The origin of these critical 
phenomena is elucidated by analyzing the structure of the phase space. An atomic cloud, containing 
the coherent modes, possesses several interesting features, such as interference patterns, interference 
current, spin squeezing, and massive entanglement. The developed theory suggests a generalization 
of resonant effects in optics to nonlinear systems of Bose-condensed atoms. 
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I. INTRODUCTION 



Low-temperature properties of dilute Bose-Einstein condensate are well described by the Gross-Pitaevskii equation 

(sec reviews [1 3]). The mathematical structmc of the latter equation is that of the nonlinear Schrodinger equation. 
This equation in the presence of a confining potential possesses a discrete spectrum of stationary states. The ground- 
state solution corresponds to the equilibrium Bose condensate. If it could be possible to macroscopically populate a 
nonground-state energy level, this would correspond to the creation of a nonground-state Bose-Einstein condensate. 
The feasibility of creating such nonground-state condensates was advanced in Ref. [4] . The properties of these excited 
self-consistent states and different ways of their formation have been previously discussed in [4-12]. A nonlinear dipole 
mode in a two-component condensate was excited in experiment [13] . 

A few words are in order to clarify the meaning of the self-consistent excited states. First of all, they should not 
be confused with the elementary collective excitations. The latter correspond to small oscillations around a given 
trap state, say around the ground state, and are described by the linear Bogolubov-De Gennes equations. While 
the self-consistent trap states are described by the nonlinear Gross-Pitaevskii equation. This is why the stationary 
solutions to the latter equation arc called the nonlinear modes. The Bose-condensed atoms possess a high degree of 
coherence [1-3] , and the related nonlinear Schrodinger equation is known to describe a coherent wave function, or the 
coherent field [14]. Because of this, the nonlinear trap modes are termed the coherent modes. The atomic densities, 
corresponding to different nonlinear coherent modes, possess rather different spatial distribution, because of which 
such modes are sometimes named as topological [3,13]. Vortices are an example of these modes. The time-dependent 
nonlinear Schrodinger equation supports the existence of solutions called solitons [15,16]. The latter are certainly 
nonlinear by their origin, but they are not necessarily identical to the nonlinear coherent modes. Only those solitons 
can be identified with such modes, which are the solutions to the stationary eigenvalue problem. 

In the present paper, we study the excitation of the nonlinear coherent modes by means of the resonance pumping 
technique suggested in Ref. [4]. To some extent, this is analogous to the problem of periodically modulated nonlinear 
optical fibers [17], as opposed to the adiabatic perturbations of nonlinear soliton equations [18]. Our aim here is to 
develop a thorough theory of the nonlinear coherent modes, analyzing their features that have not been considered in 
our previous publications [4,7]. The main new points are as follows: 

(1) A more general pumping field is considered, which makes it easier to create vortices and, what is more important, 
allows the variation of the initial phase difference related to the resonantly connected fractional populations. Varying 
this phase difference makes the dynamical picture much richer, providing several novel regimes of motion. 

(2) Necessary and sufficient conditions for the feasibility of the resonant excitation of the nonlinear modes are 
derived. These conditions impose limitations on the admissible strength of atomic interactions and on the number 

of trapped atoms. The found restriction on the number of atoms is similar to the critical number arising in the 
analysis of the stability of condensates with a negative scattering length. The ways of optimizing the found limits are 
discussed. 

(3) Temporal dynamics of fractional populations and of the related phase difference is studied in detail for all 
admissible initial conditions and for varying system parameters. Stability analysis is accomplished. Investigating the 
phase portrait determines possible dynamic phase transitions. 

(4) A special attention is paid to dynamic critical phenomena, discovered earlier [19,20] by means of numerical 
calculations, but which have not received a comprehensive explanation. The origin of these critical phenomena is 
elucidated as being due to the saddle separatrix crossing by the starting point of a trajectory in phase space. A set 
of critical lines on the parametric manifold is found. 

(5) The influence of power broadening on the long lasting resonant pumping is considered. This effect, despite well 
preserved resonance conditions on the excitation frequency, may lead to the accumulation of perturbations, resulting 
in the appearance of attenuation, limiting the admissible time of pumping. This limiting time is found to be of order 
of the trapping lifetime caused by depolarizing atomic collisions. 

(6) Interference pattern and interference current are found to display two-scale oscillations, slow and fast. These 

interference effects can be experimentally observed. 

(7) Atomic spin squeezing in the presence of multiparticle entanglement is studied. The dispersion of the population- 
difference operator is always squeezed with respect to that of the transition-dipole operator. This means that the 
fractional populations can be measured with a precision that is higher than the accuracy of measuring transitional 
characteristics, such as the relative phase. The connection of these results with the notions of state coherence and 
transition coherence is discussed. 
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II. COHERENT MODES 



The Gross-Pitaevskii equation is usually treated as an approximate mean-field equation for the broken symmetry 
order parameter [1,2]. This interpretation, however, is not unique. Here we show that this equation can be considered 
as an exact equation for a coherent wave function. The latter interpretation sheds some light on the intimate relation 
between Bose-Einstein condensate and coherence. At the same time, this allows us to give a rigorous mathematical 
definition of coherent modes which will be considered in the following sections. To suggest an accurate definition of 
the nonlinear coherent modes is of principal importance, since these are the basic object studied in the present paper. 

In the field approach, dealing with the field operators il){r,t), the standard multiatomic Hamiltonian reads 



2mo 



+ U{v,t) 



ilj{r,t) dr+ 



+ 1 j ^\r,t)^\v',t)^{v-v')i,{v',t)^{v,t)dvdv' , (1) 

where mo is atomic mass; C/(r,t), external potential; and $(r), interaction potential. The field operators satisfy 
the Bose commutation relations. Trapped atoms constitute a dilute gas, for which the interaction potential can be 
presented in the Fermi form 



$(r) = A,5{r) , A, = Anh^ ^ , (2) 

mo 

in which is an s-wave scattering length. 

Coherent states [21] are defined as the eigenstates of the annihilation operator 

V(r,t) |?7) =v{r,t) \v) ■ (3) 

Here \r]) is a coherent state and the function r]{r,t), playing the role of an eigenvalue, is called the coherent wave 
function, or, for brevity, the coherent field. The Heisenberg equation of motion for the field operator yields 

ih-^v,i)=H[^]ij{v,t) , (4) 



with the nonlinear Hamiltonian 



H[^] = - ^ + U{r,t) + Asi;Hr,t) V(r,i) . (5) 
zmo 



Taking the average of Eq. (4) over the coherent state \r]) gives the equation 



in^^r]{v,t)=H[r]]rj{v,t) (6) 
for the coherent field, with the nonlinear Hamiltonian 

H[rj]^-^ + U{r,t)+AMr,t)f . (7) 

At zero temperature, when all N atoms are condensed into a coherent state, the coherent field r]{r,t) is normalized 
to N. It is convenient to define 



r]{r,t) = VN ip{r,t) . (8) 

The external potential 

U{T,t) = U{r) + V{r,t) (9) 

consists of a stationary trapping potential U{t) and of an additional perturbing potential V{r,t). Then Eq. (6) for 
the coherent field can be written as 
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where V = V{r, t) and 



ih ^ <p{r, t) = + V) ip{v, t) , (10) 



H[^]^-^ + Uiv)+NAMr^,t)f . (11) 



Equations (6) and (10) are exact equations for the coherent fields. Such equations in physical literature are 
termed the time-dependent Gross-Pitaevskii equations. In mathematical parlance, these are the nonlinear Schrodinger 
equations [15,16]. Being complimented by initial and boundary conditions, such equations compose the Cauchy 
problem. For nonlinear differential equations in partial derivatives, the Cauchy problem possesses a unique analytic 
solution under the conditions of the Cauchy-Kovalevskaya theorem [22] . These conditions are not fulfilled for Eqs. (6) 
or (10). Hence these equations may, in general, have several solutions. To guarantee the uniqueness of solution, it is 
necessary to impose some additional constraints motivated by the physics of the problem [15,16]. In particular, it is 
necessary to specify a functional space to which the sought solution is supposed to pertain. For the studied equations, 
it is possible to proceed as follows. 

Let us consider the eigenproblem 

H[tp„] tpn{r) = En (pn{r) , (12) 

defining the stationary solutions <^n(r) labelled by a multi-index n. In the presence of a confining potential U{r), the 
spectrum {En} is discrete. The countability of the set {(^„(r)} for a nonlinear problem is not absolutely guaranteed. 
For the Gross-Pitaevskii equation with a single-well confining potential, the eigenproblem can be treated as an 
analytical continuation of that for the associated linear Schrodinger equation resulting in the limit of vanishing 
nonlinearity [10,12]. Then the set {(pni^c)} is easily countable. But in general, e.g. in the case of multiwell confining 
potentials, there appear the stationary solutions to the Gross-Pitaevskii equation, which have no linear counterparts 
[12], because of which the countability problem becomes rather nontrivial. However, for what follows, the countability 
of the set of stationary solutions to Eq. (12) is not compulsory. The main is that these solutions be appropriately 
labelled by a multi-index. The solutions (pn{^) to the nonlinear eigenproblem (12) are the nonlinear coherent modes. 
Generally, some of the levels can be degenerate, which does not hamper the consideration. Moreover, the problem 
of degeneracy, when it becomes annoying, can be avoided by adding to the Hamiltonian H[(pn] terms lifting the 
degeneracy. Such auxilliary terms may be removed after accomplishing the required procedures. The main is that 
the set {<^„(r)} be composed of all possible different solutions. This, however, does not presuppose that all functions 
from the set are mutually orthogonal, so that the scalar product 

{ipm,V>n) = J V^ir) Vnir) dr (13) 

is not compulsory a Kroneker delta. The possible nonorthogonality of the coherent modes may be due to the fact 
that the Hamiltonian li[(p] is nonlinear. In such a case, one has 

(^H[iPm]iPm, y^nj " (y^m, H[ipn]ipnj = {Em - En){ipm-, ^n) ■ 

The left-hand side of the latter equality, for the nonlinear Hamiltonian (11), reads 

which, in general, is not zero. Therefore, the modes cpm and ip„, corresponding to E^. 9^ En, are not necessarily 
orthogonal. But the modes can, of course, be always normalized, so that {^n^fn) = 1- Now, let us organize a closed 
linear envelope over the set {<^„(r)}, which, being equipped with the scalar product (13), composes a Hilbert space 

n=ZWn{r)}. (14) 

By this construction, the set {</J„(r)} is total [21], which implies that a function from the Hilbert space (14) can be 
presented as an expansion over this set. For instance, looking for a solution of Eq. (10), pertaining to the space (14), 
we may write the mode presentation 

ip{r,t) = Y,bn{t)fn{r) . (15) 
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One should not confuse the property of totahty with that of completeness. The totality of a set {(/3„(r)} means the 
possibility of presenting a function (p e H from the Hilbert space (14) as the sum (15). While the completeness of a 
set {<^„(r)} implies the existence of the resolution of unity 

^Pn ^n{r) <(r') = (5(r- r') , 

n 

in general, involving a weighting set {pn}- More details on the delicate difference between totality and completeness 
can be found in the book [21]. 

The property of totality is weaker than that of completeness. A complete basis is certainly total, while a total 
set is not compulsory complete. In our case, the set {(/?„(r)} of coherent modes, defined by the eigenproblem (12), 
may happen to be complete. A hint to this is the observation [10] that the eigenstates of the nonlinear problem (12) 
can be considered as an analytical continuation, with respect to the nonlinearity constant Ag, of the eigenstates for 
the related linear problem, which possesses a complete basis. Moreover, the completeness of the countable set of 
eigenfunctions to the one-dimensional nonlinear Schrodinger Hamiltonian has been rigorously proved [23,24]. There 
is no such an exact proof of completeness for the three-dimensional case. But we would like to stress that nowhere 
in what follows we shall use the completeness of the set {(pn{r)} and will never invoke the related resolution of unity. 
The sole thing we shall employ is the totality of the set of coherent modes, which is equivalent to the presentation 
(15). Note that for a homogeneous gas, the modes <^„(r) reduce to plane waves [25]. 



III. RESONANT EXCITATION 



In the case of a discrete spectrum, the energy levels are separated, and the interlevel distance is characterized by 
the transition frequency 

i^mn = ^ {Em " E^) . (16) 

In order to single out a pair of levels, it is necessary to connect them by a resonant transition. For this, an external 
oscillatory field 

V{r,t) =Vi{r)cosujt + V2{r)smujt (17) 

has to be imposed on the atomic system, with a frequency tuned to the transition frequency between a couple of 
chosen levels, say the levels with the energies Ei < E2, so that the transition frequency be 

U21 = ^{E2-Ei) . (18) 

Then the resonance condition reads 

< 1 (Aw = LJ- UJ21) . (19) 

The resonant excitation of nonlinear modes in Bose-Einstein condensates is rather analogous to resonant two-level 

transitions in isolated atoms [26]. The difference is that in atoms, one considers single-electron levels, while in trapped 
Bose condensate, one deals with collective coherent states of many atoms. To separate out a pair of levels, coupled 
by a resonant transition, it is necessary that the spectrum be nonequidistant. This is the case for electronic spectra 
in atoms. And, fortunately, this is also true for the spectrum of coherent modes of trapped Bose condensates, where 
the spectrum becomes nonequidistant due to atomic interactions [3,4]. It may happen, that the transition frequencies 
(16), though being different for Em 7^ are densely packed, as e.g. is the case for highly excited states of an 
optical atom [26]. Then the resonance condition (19) has to be strengthened, so that the detuning Au be smaller 
than the transition frequencies between the energy levels neighboring the selected levels with the energies Ei and £'2. 
Such a situation, because of the uncertainty relation Au ■ At ^ 1, would lead to a long time At reqiurcd to excite 
a specific mode from a densely packed spectrum. The same difficulty may arise for Bose-condensed atomic gases, 
although for the latter one has more possibility of overcoming the problem, as compared to the case of a single atom. 
This is because a single atom possesses a fixed energy spectrum, while the spectrum of nonlinear coherent modes 
can be modified by varying the number of particles A'^ or the strength of atomic interactions, e.g., by means of 
Feshbach resonance. At the same time, these interactions render the system nonlinear, which leads to some important 
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differences in the behaviour of Bosc condensate, subject to the action of a resonant field, as compared to optical 
resonance in atoms. Comparing these two types of resonance, we may say that in a single atom, one can induce a 
linear resonance, while in a collective of coherent atoms, there occurs a nonlinear coherent resonance. 

The behaviour of the coherent field (15) is described by the nonlinear equation (10), where the alternating field (17) 
is in resonance with the transition frequency (18) of a chosen interlevel transition. Earlier, we have considered [4,7] 
only the transitions from the ground state to an excited collective state. But if one is able to transfer atoms from the 
ground state to an excited level, then by applying one more alternating field, with a different resonant frequency, it is 
possible to transfer atoms from one excited level to another higher level. Therefore, we consider the general situation 
of an arbitrary pair of collective levels coupled by a resonant transition. 

The coefRcients 6„(t) in the mode presentation (15) can be written as 

bn{t) = c„(i) exp ( - i S„ i ) , (20) 



h 

with En being the n-mode energy. In order that different modes could be separated, it is necessary that the coefficient 

Cn (t) be a slow function of time, as compared to the fast oscillations of the exponential ex.p{—-^Ent). This implies 
the validity of the necessary condition for mode separation 



h 

En 



dCn 



dt 



«: 1 . (21) 



Inequality (21) reminds us the slowly-varying amplitude approximation, so common in optics [26]. Under this condi- 
tion, we may try to get a nontrivial mode picture, following the way accepted in optical resonance. For this purpose, 
we substitute the mode form (15) into the evolution equation (10), multiply the latter by i^* (r) exp(|i5„t), integrate 
it over r, and average over time, treating c„(i) as quasi-invariants. The mathematical foundation of this procedure is 
based on the averaging technique [27]. In this way, we meet two transition amplitudes, one that is caused by atomic 
interactions 

amn = ^sj J |v'm(r)|' [2|(^„(r)|2 - |^„(r)|2] dv , (22) 
and another due to the alternating field, 

/3mn = ^J flir) [V,{r) - iV2{r)] ^„(r) dr . (23) 

Recall that in optics one has only the transition amplitTidc of an external resonant field defining a Rabi frequency. 
Here, because of atomic interactions, there appears an additional amplitude (22). Note that for different energy levels, 
the amplitude (22) may have different signs, for a fixed sign of the interaction parameter As- 

The condition for mode separation (21) is necessary, however not sufficient. For employing the averaging technique 
[27], we need that the transition amplitudes be smaller than the related transition frequencies, 

< 1 . (24) 

This mathematical requirement for the applicability of the averaging technique has a transparent physical meaning — 
all perturbations, either produced by internal interactions or by external fields, must be smaller than the corresponding 
transition frequencies, in order that the latter be well defined as such. Again, the situation here is completely analogous 
to optics, where one requires that the Rabi frequencies be smaller than the transition frequency, in order to avoid the 
power broadening spoiling the resonance picture. In our case, the difference with optics is that we have, in addition to 
the external transition amplitude (23), an internal amplitude (22). The latter appears because of the interaction term 
in the nonlinear Hamiltonian (11). This term depends on time through the hmction ip{r,t). Therefore, one could 
naively decide that as soon as the coherent field (15) changes from one mode to another, the resonance condition 
would be broken. However, this does not happen if the restrictions (24) are hold true. Then, despite the variation of 
the interaction term, the transition frequencies are always well defined, making it possible to tune a well preserved 
resonance. 

Restrictions (24) are stronger than inequality (21). The latter is only necessary, while the former are necessary and 
sufficient for mode separation. Mathematically, conditions (24) allow us, when employing the averaging technique, to 
use the integral 

1 

hm — / ex.p{tUJmnt) dt = 5mn ■ 
T^oo 1 Jq 
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As a result, we come to the set of equations 



dcn \ ^ 



2 ^n2Pl2^ie 



-iAujt 



(25) 



for the functions c„ = c„(t), satisfying the normaUzation condition 

n 

Prom Eq. (25) it immediatefy follows that, if n ^ 1, 2, then one has 



Cn{t) = c„(0)exp ■ 



E 



J' 

Jo 



\Cmit')f dt' 



(26) 



Supposing that at the initial time no levels, except n = 1, 2, are populated, that is c„(0) = for n 1, 2, we get from 
Eq. (26) that 



Hence, the normalization condition is 



c„(t) = (n^l,2;i>0). 



|ci|2 + |c2|2 = l. 



And the set of equations (25) reduces to two equations 

i^ = ai2|c2|^ci + i/3i2C2e*^"*, 
This reduction is valid provided that the conditions 



« = a2i|ci| C2 + -Pi2Cie 



iAujt 



ai2 


«1, 


0:21 


« 1, 


/3l2 


< 1 


t021 




L021 


W21 





(27) 
(28) 

(29) 
(30) 



are satisfied. 

It would be tempting to postulate from the very beginning that, under the resonance condition (19), only two 
modes are involved in the process. However such a simplified approach would leave doubts of the possibility to reduce 
the consideration to an effective two-level system. Therefore, we felt it necessary to present an accurate mathematical 
derivation of the reduced equations (29). Moreover, in the process of this derivation, we found conditions (30) defining 
the region of applicability for the resonant two-lcvcl picture. 

By acting on the trapped Bose gas with several resonant fields, we could separate out not two but several coherent 
modes, reducing the consideration to a finite-level system. The possibility of creating from Bose condensate the 
resonant finite-level systems is of great importance. This allows for a wide variety of different applications, analogous 
to those for resonant finite- level atoms in optics [26,28]. Since the properties of resonant Bosc-condcnscd systems are 
rather different from those of single resonant atoms, wo may expect the appearance of novel effects, such as dynamic 
critical phenomena [19,20]. In what follows, for brevity, the effect of resonant formation of nonlinear coherent modes 
in Bose-condensed gas will be called coherent resonance. 

Note that a resonant finite- level system of coherent modes is principally different from the two-mode models accepted 
for describing Bose condensates in stationary double- well potentials [8,29-30]. Though some equations in our case 
may resemble some expressions in the case of the double-well potentials, this similarity is rather formal, being due to 
the simple fact that all finite-level systems share some common properties in their mathematical structure. However 
the physics of a system, subject to the action of a resonant alternating field, drastically differs from what happens in 
stationary double wells. 



IV. TRANSITION AMPLITUDES 



The set of inequalities (19) and (30) provides us with the necessary and suflScient conditions for realizing the coherent 
resonance. The frequency and the amplitude of the alternating external field (17) can always be chosen in order to 
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satisfy the resonance condition (19) and the last of the inequalities (30). Hence the main concern is connected with 
the validity of the first two inequalities (30) for the internal transition amplitudes (22) due to atomic interactions. To 
understand better what is the restriction imposed by these inequalities, let us consider a harmonic cylindrical trap 
with the radial frequency LOr = (^x = and axial frequency w^, denoting the aspect ratio 



(31) 



Introduce the dimensionless cylindrical variables 



' X ~ ' y r 



Z=f Vr = \ , (32) 

measured in units of the oscillator length ly. In the dimensionless coupling parameter 

= 47r ^ iV , (33) 

as is an s-wave scattering length and N is the number of particles. The nonlinear Hamiltonian, measured in units of 
Wr, is 

U = -)^^ + \{T'' + v^z^)^a\^\\ (34) 

where 

V'(r,^,z) = ^^/^(r) 

is a dimensionless wave function and 

V2--^ - — — — — 

In this notation, the eigenproblem (12) reads 

(r, (p, z) , (35) 

where n = 0,1,2,... is the radial quantum number, m = 0, ±1,±2, ... is the azimuthal quantum number, and 
j = 0, 1, 2, ... is the axial quantum number. 

The eigenproblem (35) can be solved by employing the optimized perturbation theory [33-35]. For this purpose, 
we take as an initial approximation the Hamiltonian 

^o = -^V^ + i(nV+«V) , (36) 

containing two trial parameters, u and v, playing the role of effective dimensionless frequencies. This Hamiltonian 
possesses the eigenfunctions 



'"'"""'^'^''%Hexp(-|r^)iH (^,2^ x 



{n+\m\y. 



X —= - , exp - - 2 HA^ z) , 

^/2^ W/ ^/277! V 2 ) ^^"^ " 

where L™(-) is a Laguerre polynomial and ^/j(-) is a Hermit polynomial. Then, invoking the Rayleigh-Schrodinger 
perturbation theory, we may find the higher approximations. Thus, for the spectrum, we get the sequence {Ek{g, u, v)} 
of approximations, with A; = 0, 1, 2, . . . being an approximation order. For instance, in the first order, 

Ei{g,u,v) = ^(u+^\ +^(v+'^\ + uV^ ginmj , (37) 
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where the dependence of the energy on the quantum numbers is not explicitly marked in order to avoid cumbersome 
notation; the quantum numbers enter through the combinations 

p = 2n+\m\ + l , q = 2j + l 

and through the integral 

Inmj = I \tpnmj{r, ^fi, z)f rdr dip dz . 

Then the trial parameters u and v are to be transformed to control functions U]^{g) and Wfe(g), which play here the 
role of control frequencies, such that the sequence {E]^{g,Uk{g),Vk{g)y} be convergent [33]. This transformation is 
accomplished by means of an optimization condition. Limiting ourselves by the first-order approximation, we may 
write the optimization condition as 

'^^TT + ^^Tt) M9.u,v) = Q. (38) 



du dv 

In view of expression (37), this gives the equations 



p|i-i) + ^,/^ = o, ,fi-Il!Vr:^ = o (39) 



pv y q \ J pvy/vq 

for the control frequencies u = u{g) and v = v{g), where the notation 

G = 2py/q Inmj gi^ (40) 
is used. Substituting these control functions into Eq. (37), we obtain the optimized approximant 

E(g) = Ei(g,u(g),v{g)) (41) 

for the energy spectrum. Recall that the energy (41) as well as the control frequencies u = Umnj and v = Vnmj depend 
on quantum numbers. 

The spectrum (41), for arbitrary values of g, is defined by the system of equations (37) and (39). In the limits 
of weak coupling and strong coupling, we can derive explicit asymptotic expansions. Thus, for weak coupling, when 
5 — > and G ^ 0, the control frequencies are 

G Gu 



2p\qvY/^ ' ~ 2p{qvf/^ ' 

And the energy spectrum is 

Ec±ao + aiG (G ^ 0) , (42) 



where 



2 ' ' 2p{qu)y^ 
In the limit of strong coupling, the control frequencies are 

p qv"^ 



Q2/5 ' - Q2/5 ■ 

And for the energies, we find 

E ~ boG^^^ + 6iG-2/^ (G ^ oo) , (43) 

where 

5 2p^ + {qi^f 
^0 = 4' ^1 = 4 • 
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The difference l)etwcen two selected energy levels defines the transition frequency (18). 

If we consider transitions between the ground-state level {n = m = j = 0) and an excited level, then calculating 
the internal transition amplitudes (22), we meet the integral 



'P, z)\'^\tpnmj{r, (p,z)\'^ rdr dip dz , 



(44) 



in which uq = 'Uooo(fl') and t;o = t^ooo(fl')- For the transition amplitudes (22), corresponding to transitions between the 
ground state and an excited mode, labelled by the quantum numbers n, m, and j, we obtain 



'Ur. 



^ 71711 J 



where Iq = /qoo- The first two inequalities in Eq. (30) for the internal transition amplitudes now read 



nmj 



< 1 



< 1 . 



Introduce the quantities 



(45) 



(46) 



Oir, 



3,0 ) 



:(ao,r; 



,,o) 



(47) 



Here 



is an average transition amplitude, and Snmj plays the role of an internal detuning induced by atomic 



interactions. By direct numerical calculations of Eqs. (45) and (47), we have checked that \Snmj / <^nmj\ ^ 0.1. In 
appendix, we demonstrate calculations for several first excited modes. 

The analysis of inequalities (46) shows that their validity imposes a constraint on the renormalized coupling (40), 
requiring that the latter be outside the region of convergence of the strong-coupling expansion (43). Explicitly, this 
yields the condition 



[2p2 + jqi^)^]'/^ 

Inmj 



(48) 



The constraint (48) on the coupling parameter (33) implies that the number of particles in the trap has to be less 
than the limiting number 



[2p^ + jqu)^]^/^ 
56nup^ Inmj 



(49) 



The value of Nq depends on the quantum numbers, thus, being different for different modes. For the ground state, 
Eq. (48) becomes 

\gu\<{2 + u'f' , 
and, respectively, for the limiting number (49), we get 

The limiting number increases for the higher excited modes. Thus, it grows as 

TVo - (2n + \m\ + 1)^/2 , iVo - (2j + 1)^ 

according to whether the radial or axial modes are excited. For highly excited modes, one could invoke an optimized 
expansion in powers of h [36] . The number A^o is also very sensitive to the trap shape, depending on the aspect ratio 
(31). For the disk-shape (i^ 3> 1) and cigar-shape (i^ <C 1) traps, Nq is larger than for a spherical {v = 1) trap. 

It is interesting that the limiting number of particles for an excited coherent mode is close to the critical number 
of atoms for the stability of a trapped Bose gas with attractive interactions. Hence, such gases can also support 
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coherent modes. Examples of atoms with a negative scattering length are ^Li (see review [37]) and ^^Rb (see Ref. 
[38]). The recent developments in utilizing Fcshbach resonances have shown the ability to tune the scattering length 
to almost any desired value, including the ability of changing its sign [3], which makes it possible to add to these 
two species many other atoms. The critical number of trapped Bose-condensed atoms with attractive interactions 
has been estimated earlier [4,39-41]. Such trapped atoms form a metastable state, which decays because of quantum 
tunneling, which however, is very slow [42 44], if the number of atoms is less than the critical munber. There was a 
proposal [45] for stabilizing Bose condensate with attractive interactions by driving a quadrupole collective excitation. 
Our consideration above suggests that it could, probably, be also possible to stabilize such condensates by transferring 
them in an excited coherent mode with the help of the coherent resonance. 

To estimate the limiting number of particles in a coherent mode, consider a cigar-shape trap, as was used [46] 
for condensing ^^Na. With the radial frequency Wr = 1527 Hz and axial frequency = 11 Hz, the aspect ratio 
is v = 0.007. The scattering length of ^^Na is a,, = 4.498 x 10~^ cm. Since the oscillator length in this case is 
Ir = 1.345 X lO^"' cm, then Ir/ag = 300. The coupling parameter has to be <? < 340. Therefore, the limiting number 
A^o ~ lO''- If one would take a very long cigar-shape trap with the aspect ratio v = 0.001, then g < 2400 and the 
limiting number could be as large as A^o ^ 10^. 



V. PHASE PORTRAIT 



To study the dynamical behaviour of the system under coherent resonance, it is convenient to use the population 
difference 

s = \c2f - \cif , (51) 
which varies in the interval — 1 < s < 1. Then, in view of the normalization condition (28), we have 

112 1 — S I |2 1 + S 

jCll \C2\ 

Taking this into account, we may write 



ci = Y exp |i ^TTi -h ^ f j I , C2 = y' exp |i (w2 - ^ * j } ' (^2) 

where tti = Tri{t) and tt2 = T'2{t) are real- valued phases. 

To simplify the notation, we introduce the average interaction amplitude 

a = ^ (ai2 + a2i) , (53) 

which is real. Noticing that the transition amplitude (23), due to the resonant field, is, in general, complex- valued, 
we can present it as 

P,2=pe'\ P=\Pi2\. (54) 

Also, we define the effective detuning 

6 = Au}+^ (ai2 - Q!2i) . (55) 
And, finally, we shall need the phase difference 

X = 7r2 - TTi -I- 7 , (56) 

whose initial value Xq can be varied by changing the spatial dependence of the resonant field. With this notation, 
Eqs. (29) transform to the system of equations for the population difference 

^ = -p Vl-s^ sin a; (57) 
at 

and for the phase difference 
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— = as-\ — , cos a; + d . (58) 

dt y/l - s2 



These are the equations that can be written in the Hamiltonian form 

ds _ dH dx _ dH 
dt dx dt ds 

with the Hamiltonian 



H{s,x) = - -/? - s2 COSX + 5s . (59) 

Let us consider the stationary sohitions to Eqs. (57) and (58) in the rectangle of the variables — 1 < s < 1 and 
< a; < 27r. These solutions arc given by the equations 

s^ + 2es^ - {l-l? - e^) s^ -2£s- = Q , sina; = 0, (60) 

where the notation 

6^^, e^^- (61) 
a a 

is used. In what follows, for simplifying formulas, the dimensionless detuning e will be treated as a small parameter, 
|£| <C 1. As we have checked by direct numerical calculations, the quantity (ai2 — a2i)/a is small, and in addition, it 
can always be conipcmsatcd by Aw, so that 6 /a be as small as required. 

The phase-portrait picture depends on the value of the parameter b characterizing the amplitude of the resonant 
field. This parameter, according to notation (61), can be either positive or negative, depending on the sign of the 
interaction amplitude a. The latter can be negative in the case of attractive interactions. Also, it can become negative, 
as follows from definition (22), even for repulsive interactions, when the resonantly connected levels are such that 

2(bm|^ Wn\'^) < (bmr> ■ 

We calculated numerically the scalar products entering the above inequality and found that, when exciting the lower 
modes, the amplitude a does not change its sign. For this to occur, the resonantly connected modes must be 
energetically strongly separated. 

When 6^ > 1, there are three stationary solutions to Eqs. (57) and (58), 

= I = S3 , xl=Q, xl = 2-K , 

S2 = - f , xl=T,. (62) 
For < 6 < 1, in addition to the stationary points (62), there appear two other fixed points, 

S4 = + z — , a;^ = TT 



1-62 



S^ = -VT^+^— ^, xI=t:. (63) 

The points (62) and (63) are the solutions to Eq. (60), which are simplified taking account of small detuning |£| <C 1. 
When — 1 < 6 < 0, the phase portrait contains seven fixed points, those given by Eq. (62) plus 

54 = Sg , X4 = , Xq = 27r , 

sl = s*r, x*^=0, x*r = 2n, (64) 

where s| and S5 are the same as in Eq. (63). A qualitative change of a phase portrait is what in dynamical theory 

called a dynamical phase transition. In our case, as follows from Eqs. (62) to (64), there are such dynamical phase 
transitions at the values of 6 = and = 1. One may notice that the phase portraits for 6 > and 6 < are similar 
to each other, just being shifted with respect to x by n. 
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The Jacobian matrix associated with Eqs. (57) and (58) is composed of the elements 



Ji2 = -b^/l 



b cos X bs 

•^21 = 1 + (1 _ g2)3/2 ' -^22 = - ^^^==f Smx . 



The eigenvalues of this matrix are given by the equation 



.2 



= 5- (s^ sin^ X — cos^ x) —b y/l — s"^ cos a; . (65) 

1 — S'* 

These eigenvalues are to be evaluated at the fixed points listed above. Then for 6^ > 1, we have 



J^t = ±i VHb+l) = , Ja = ±* Vb{b-l) , (66) 
which shows that all fixed points (62) are the centers. When < 5 < 1, then 



Jt = ±i Vb{l + b) = Jf , = ± - b) , Jf = ±i = Jf , (67) 

from where if follows that the fixed point (s2,a;2) is a saddle, while all other points in Eqs. (62) and (63) are the 
centers. For — 1 < 6 < 0, the fixed points (64) are similar to those (63), being shifted by tt with respect to x. All 
points remain the centers, except (sj, .Tj) which is a saddle. 

The trajectory, passing through a saddle point, is the saddle separatrix, wliic;li separates the phase plane onto 
qualitatively different regions of motion. In the considered case, the separatrix is given by the equality 

H{s,x)=H{sl,xl), (68) 

which, in the approximation linear in e, results in the separatrix equation 

„2 



2 — 6 \/l — cos x + es = b . (69) 

This equation defines two separatrices the upper and lower ones. If the initial point (sq, .tq) of a trajectory lies below 
the lower separatrix, then the variation of s{t) is always limited from above. In the same way, if the starting point 
of a trajectory is above the upper separatrix, then the change of s{t) is always limited from below by the separatrix. 
Since varying the parameter b moves the separatrices, which is clear from the separatrix equation (69), the whole 
phase picture essentially changes. The transformation of the phase portrait with varying the parameter 5, related 
to the increase of the transition amplitude due to the resonant field, is illustrated in Fig. 1, where only the most 
characteristic pictures are presented. The detuning e in these figures is set zero. A finite detuning £ ^ makes the 
pictures slightly asymetric with respect to the line s = 0, but the overall phase portraits are similar, because of which 
we do not present them here. 

The structure of the phase portraits shows that, at small b, there are a number of trajectories always lying in the 
upper part of the plane if their initial points were above the upper separatrix. In particular, if the upper coherent 
mode was initially highly populated, and the applied resonant field is weak, the mode will retain its high population. 
In the extreme case, if at the initial time all atoms are in the upper coherent mode, and there is no external field, so 
that So = 1 and 6 = 0, then s{t) = 1 meaning that atoms remain in the upper mode. Such a setup can be achieved in 
two ways. For instance, by transferring atoms to the upper level through an intermediate third higher mode. Or, even 
simpler, just by applying, first, a resonant field sufficient for transferring atoms from the ground state to the chosen 
coherent mode and then switching off this field. Thus atoms remain locked in this upper mode. This mode-locking 
effect was described in Ref . [4] . Mathematically, it is analogous to the self-trapping of atoms in one of the wells of 
a stationary double- well potential [30-32,47]. However, the physics of the coherent resonance occurring under the 
action of a resonant external field is absolutely different. Moreover, the mode-locking effect is not the most interesting 
thing happening in the process of the resonant excitation of coherent modes. 

Varying the amplitude of the resonant field, that is, the parameter b, in the dynamic behaviour of fractional 
populations there appear dynamical critical phenomena [19,20], arising on a critical line of the parametric manifold 
{b, e}. The qualitative behaviour of fractional populations drastically changes when the critical line is crossed. As was 
shown [19,20], for a stationary system, obtained by time-averaging the considered dynamical system, there appear 
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critical phenomena, typical of phase transitions and occurring on the same critical line, where the dynamical critical 
phenomena are observed. Thus, we may define a pumping capacity, playing the role of an effective specific heat, and 
an effective susceptibility, which diverge on the critical line [19,20]. The origin of these dynamical critical phenomena 
has not been completely understood. But now, invoking the analysis of this section, we clearly see what happens. 

Suppose, we consider a trajectory starting at the point {so,a;o}. Strengthening the applied resonant field means 
increasing the parameter b. But changing b implies moving the separatrix given by Eq. (69). When the separatrix 
crosses the initial point {so,xq}, the trajectory passes from one region of the phase plane to a qualitatively different 
region, because of which the dynamical behaviour changes qualitatively. Since, in general, we have two independent 
parameters, the amplitude b and detuning e, the effect of separatrix crossing can be due to the variation of any of 
these parameters. A tiny change of one of these can shift the trajectory to a different part of the phase plane. The 
condition for the separatrix crossing defines, on the parametric manifold {b,s}, the critical line 



Y -be yl- si cosxo+£cSo = \bc\ ■ (70) 

Actually, for each set of initial conditions {sq,xo} there exists its own critical line defined by Eq. (70). For the variety 
of initial conditions, there is an infinite bunch of critical lines. 

For example, let us consider the initial condition sq = —1, Xq = 0, when all atoms at the initial time are in the 
ground mode. Then the critical line (70) reduces to 

l^cl +£c = ^ ■ 

In the case of zero detuning, this defines two points b,, = ±0.5. Assume that our concern is the behaviour of atoms with 
repulsive interactions, such that 6 > 0. And let us set £ = 0. The population difference x{t) exhibits dramatic changes 
when b crosses the critical point be = 0.5. This is illustrated in Fig. 2, where the change of b around the critical point 
be occurs in the fifth decimal digit, so that the variation of b from Fig. 2(a) to Fig. 2(b) makes Ab = 0.00001. Slightly 
below the critical point be, the population difference s(t) oscillates between —1 and 0, while the phase difference x{t) 
monotonely diminishes. Just above the critical point, the oscillation period of s{t) doubles, as well as the amplitude 
of s{t), as now s{t) varies between —1 and +1. The form of s{t) also changes acquiring a specific cusp in the middle 
of the period. The phase difference x{t), instead of being monotonely diminishing, becomes oscillating. 

If we are interested in the behaviour of atoms with 6 < 0, say with attractive interactions, then similar dramatic 
changes happen at the critical point be = —0.5. This is shown in Fig. 3. Again, the change of b around the critical 
point be occurs in the fifth decimal digit, as for Fig. 2. After crossing the critical point be, the period and amplitude 
of s{t) double, and the middle-period cusps appear. The shapes of both s{t) and x{t) essentially change. Note that 
although these dynamical critical phenomena for 6 > and 6 < are similar, there is a difference in the behaviour of 
x{t). For 6 < 0, the phase x{t) is an oscillating function below as well as above the critical point. While for 6 > 0, 
the phase difference x{t) below the critical point is a monotonely diminishing function. 

In Figs. 2 and 3, time is dimensionless (measured in units of a~^). The return to dimensional time is accomplished 
by the substitution t ^ t/a. A small detuning e ^ slightly curves the lines but does not make qualitative change. 



VI. POWER BROADENING 



The feasibility of coherent resonance imposes a limitation on the value of the internal transition amplitude, that 
is, as discussed in Sec. IV, on the admissible number of particles. In addition to the amplitude restriction, there is 
one more limitation having to do with temporal power broadening. Even if the alternating external field is well tuned 
to the resonance with a chosen transition frequency, nevertheless, there exist nonresonant transitions to other levels. 
Though the probability of nonresonant transitions, at each given instant of time, is very small, their influence on the 
process of resonant excitation can increase with time, becoming essential for long time intervals. To estimate this 
influence, we return to the evolution equations (29). 

It is convenient to introduce the quantity 

h = 2ctc2e'(^'"*+T) = e'^ , (71) 

playing the role of a ladder variable in optics [26,28]. Here s = s{t) and x = x{t) are the same as in Eqs. (51) and 
(56). Define the collective frequency 

n = as + S . (72) 
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Here and in what follows, wc use the notation for a, (3, and 5 as given in Eqs. (53), (54), and (55), respectively. 
Taking account of nonresonant transitions can be done by including in the evolution equations an inhomogeneous 
broadening modelled by a random variable [28]. Then, Eqs. (29) can be rewritten in the form of the stochastic 
equations 

^=i{n + ^)h + il3s , (73) 



(74) 

where ^ = ^(i) is a stochastic variable modelling nonresonant transitions. 

In analysing Eqs. (73) and (74), we assume, for simplicity, that the transition amplitude due to the external field 
is weaker than that caused by internal atomic interactions. 



< 1 . (75) 



The stochastic variable can be presented as Gaussian white noise [48] , with the stochastic averages 

« m > = , < i{t) i{t') > = 273 5{t - t') . (76) 
The noise width can be approximated as 

73 « — , (77) 

SO that, according to inequalities (30), 73 <C u;2i- 

Equations (73) and (74) can be solved by employing the scale separation approach [49-51], which is a generalization 
of the averaging technique [27,52] to stochastic differential equations. Because of inequality (75), the fimction h = h(t) 
can be treated as fast, compared to the slow function s = s{t). The latter is a temporal quasi-invariant with respect 
to h. In this case, the solution to Eq. (73) reads 

h = /loe'^* exp |i ^{t') di'l + il3s e*"(*-*') exp |i £ ^{t") dt"^ dt' , (78) 

where ho = h{0). Averaging Eq. (78) over the stochastic variable, with the use of Eqs. (76), we have 

« /I » = ^0 + -^"j e(^"-^3)* _ . (79) 

Substituting the solution (78) for the fast function into Eq. (74) for the slow function and averaging the right-hand 
side of Eq. (74) over the stochastic variable and over time, we obtain the equation 

ds__ fi^ss 
dt (as + 5)2+71 

for the guiding center of the slow function. An approximate solution to this equation can be written as 



«oexp ( - ^2 , o ) • (81) 



This shows that the population difference (81) attenuates, with the characteristic time 

_ {asp + 6f + 7I 

= — WTs • 

Prom here it follows that the process of resonant pumping is not essentially influenced by power broadening only 
during times smaller than tc- After the time (82), the two-level resonant picture will not be adequate, since the 
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neighboring nonresonant levels will be essentially involved in the process. The characteristic time (82), in view of Eq. 
(77), is 

Taking account of the inequalities \6/a\ <C 1 and \P/a\ <C 1, and setting |so| ~ 1, simplifies the time (83) to 

*c-^. (84) 



Finally, under conditions (46), Eq. (84) reduces to 



ic ^ ^ • (85) 



To estimate tc we may take /3 ~ 0.1a and a ~ 0.1a;2i, with a transition frequency, typical of magnetic traps [3], of 
order W21 ^ 10^ — 10^ Hz. This gives tc 10 — 100 s, which is a rather long time, comparable or longer than the 
lifetime of atoms inside a trap [2,3]. 

In this way, the temporal limitation resulting from the power broadening does not impose too severe restrictions 
on the procedure of the resonant excitation of coherent modes. Moreover, as is discussed in Sec. V, we do not need 
to pump too long but we may stop pumping as soon as the desired mode is populated, which happens during the 
time of order 27r/0. The latter is of order 0.1 s. Hence, the time necessary for populating a chosen coherent mode is 
essentially less than the critical time tc, after which power broadening would spoil the resonant picture. 

VII. INTERFERENCE PATTERNS 

In the process of the resonant excitation of coherent modes, only two modes are involved in the dynamical picture, 
while the population of all other nonresonant modes remains negligibly small. This allows us to present the system 
coherent function as the sum 

ip{r,t) = ipi{^,t)+M^,t) (86) 

of the terms 

fi{r,t) = c,(t)(^,(r)exp (- E,t^ (87) 

corresponding to the related modes i = 1,2. From here it follows that there should exist a spatial interference between 
the modes, similarly to the interference of different atomic components in a binary mixture of two Bose condensates 
[53]. In our case, the interference becomes possible because different coherent modes possess qualitatively different 
spatial shapes. Thus, we may define [3] the interference pattern 

Pint {r,t) = p{r,t) - pi{r,t) - p2{r,t) , (88) 

in which 

p{r,t) = \ip{r,t)W p^ir,t)=\ip,{r,t)\\ (89) 
With the use of the dimensionless real function 

i;i{r,z) = l',/^^i{r)e-'^^^ , (90) 
in which Ij- is the oscillator length and rrii is a winding number, the pattern (88) reads 

where s = s{t) is defined by Eq. (51) and the notation 



P^nt{r,t) = ^ COS^ , (91) 
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$ = (m2 -mi)<fi + x{t) - 7 - {u)2i + Aco)t 



(92) 



is introduced. This interference pattern can be experimentally observed either inside the trap by moans of light 
scattering or by freeing atoms from the trap and observing their free evolution, as is discussed by Sinatra and Castin 
[53]. 

The existence of two modes in a trap, with different spatial shapes of the related wave functions, leads to the 
appearance of atomic current inside the trap. This effect is analogous to the Josephson effect, which is usually 
considered for two atomic clouds sitting in a double-well potential, so that the clouds are separated by a potential 
barrier [29- 32]. However, as is suggested by Leggctt [54]. Josephson oscillations can exist between two interpenetrating 
populations, not separated by any barrier. In this case, one calls it the internal Josephson effect [55]. Such a tunneling, 
involving no potential barriers, is also called quantum dynamical tunneling [56,57], which is actually just a current 
between two modes representing bound states. 

The interference current is 



jmt(r,i) = j(r,t) - ji(r,f) - j2(r,i) 



(93) 



where 



n 



n 



j(r,i) = — Im(/p*(r,i)VV(r,t) , ji(r,t) = — Im v?*(r, i)V(pi(r, i) 



mo mo 
with i = 1,2. Employing the notation in Eqs. (90) and (92), we find 



jmt(r, i) = 



?tV'lV'2 

2mo/4 



d_ 

^ dr 



sin $ + 6;, 



mi + m2 



cos$ 



(94) 



(95) 



Both the interference patterns (91) and interference current (95) experience temporal oscillations on two time 
scales. One corresponding to the transition frequency u;2i) which is typical of the Josephson effect. And another scale 
is related to the temporal variation of s{t) and x{t), which is conneted with the transition amplitudes a and /?. Since 
a,/3 -C a;2i, the Josephson oscillations are fast, as compared to the slow change of s{t) and x{t). The fast Josephson 
oscillations, being modulated by the slow variation of s{t) and x{t), yield the collapse-revival picture typical of some 
two-level systems [26]. Here the modulation is mainly due to the temporal variation of s{t). If one stops applying the 
resonant field, so that /3 = 0, then s{t) = const, and the slow modulation disappears. 



VIII. SPIN SQUEEZING 

The considered atomic system is essentially nonlinear. Therefore, one might expect that some kind of squeezing 
effects could arise. Similarly to squeezed states of light [58,59], one can introduce squeezed atomic states [58]. This is 
usually considered for the case of two-level atoms, each of which possesses two internal states. Finite level systems, 
as is known, are conveniently described by means of spin operators, because of which atomic squeezing is commonly 
called spin squeezing. In order to emphasize that the spin operators, employed for describing finite-level atoms, are 
not actually the operators representing real spins, but rather are convenient mathematical tools, one also uses the 
terms of dipole squeezing [58,60] or pseudospin squeezing [61]. In general, one may define squeezing for other operators 
from a Lie algebra [61,62]. Atomic squeezing is directly related to the radiation field squeezing [63]. And vice versa, 
squeezed atomic states can be created by irradiating atoms with squeezed light [64 67] or with light combined with an 
alternating magnetic field, as in the process of continuos quantum measurement [68] . Generally, atomic squeezing can 
be achieved for both bosons as well as for fermions [69,70]. Squeezing in Bose-Einstein condensates was considered for 
two-component mixtures [71], for atoms with two internal states [72], which is equivalent to a two-component mixture, 
and for atoms in linked mesoscopic traps formed by an optical lattice [73], which is equivalent to a multicomponent 
mixture. There is a real potential for several practical applications of squeezed atoms, e.g., for atomic spectroscopy 
and atomic clocks [74], for atom interferometers [75], and, probably, for quantum computation [76]. 

To consider spin squeezing in our case, we may notice that the transition dipole (71) and population difference (51) 
can be presented as the statistical averages 

h=^ <S_>, s=^ <S,> (96) 

of the collective spin operators 



17 



N 

S^ = J2sr, S± = S,±iSy, (97) 

i=l 

where a = x,y,z and 5f corresponds to a 1/2-spin operator. The evolution equations (73) and (74) can be obtained 
by averaging the Heisenberg equations of motion for the spin operators (97), with the effective Hamiltonian 

H,ff = ^iS.+S+)-SS.- ^Sl. (98) 

Atomic squeezing is defined in a way similar to the squeezing of light, being based on the quantum-mechanical 
uncertainty relations. For any two operators A and B, not necessarily Hermitian, the Heisenberg uncertainty relation 
reads 

A2(^) A2(B) >\\< [A, B] >|^ , (99) 

where the dispersion is 

A^{A) =< A+A > - \<A>f. (100) 
One says that A is squeezed with respect to B if 

A^iA) <1\<[A,B]>\. 

This suggests to introduce the squeezing factor 

Now, for A and B, we may take any spin operators. However, taking separately S^, with a = x,y, z, may involve the 
so-called trivial squeezing due to rotation [77]. To avoid this, we consider the operators Sz and S±, for which we have 

A\Sz) A^{S±) > J I < 5± > p . (102) 

Squeezing of with respect to S±, is defined by the squeezing factor Qs^s± = Qz, which in view of definition (101), 
is 

The squeezing factor (103), taking into account that \ < S± > =< Sx >^ + < Sy >^, becomes 

V< 5X >2 + < Sy >2 ^ ' 

which has the form used in Refs. [72,74,77]. One can say that Sz is squeezed with respect to S± if A'^{Sz) < |] < 
S'i > ], that is, Qz < 1. Calculating the expressions 



S2 



we obtain the squeezing factor 



Qz = Vl^, (105) 



where s = s{t) is the population difference satisfying the evolution equations discussed above. As is seen, since 
< < 1, the factor (105) is almost always less than one, except for s = 0. Hence Sz is squeezed with respect to 
S±. In other words, the dispersion A'^{Sz) is almost always smaller than the dispersion A'^{S±). This means that the 
population difference and, respectively, the fractional populations can be measured with a better accuracy than the 
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transition dipole or the relative phase between population amplitudes. The characteristic temporal behaviour of the 
squeezing factor (105) is shown in Fig. 4, with s(t) found numerically from the evolution equations. 

It is worth noting that the treatment of atomic squeezing by means of spin operators is often used in the quantum 
optics language [58-61]. This is equivalent to the formation of squeezed states generated by the Bogolubov canonical 
transformation for quasiparticles in multicomponent (spinor) condensate. These two pictures are equivalent mathe- 
matically and are caused by the same physical reason, by the existence of nonlinear interactions in a multicomponent 
system. 

In the theory of nuclear magnetic resonance [78,79], one distinguishes the state coherence, when s = ±1, and the 
transition coherence, when \h\ = 1. This terminology can also be applied to our case. Since here wc have = 1 — s^, 
then the state and transition coherences are complimentary to each other. And if Qz < 1, this means that the state 
coherence can be better controlled than the transition coherence. 

Atomic squeezing is usually neighbours with atomic entanglement [72,80,81]. Entanglement is simply the 
Schrodinger's name for superposition in a multiparticle system. In the system of trapped atoms or ions with two 
or more internal states, multiparticle entanglement can appear [81]. This concerns as well Bose-Einstein condensates 
composed of atoms having internal states [72,82]. The Bose-condensed atomic cloud, subject to the coherent resonance 
we study here, is also in a multiparticle entangled state. Really, the multiparticle density matrix, for the case studied, 
is 

N N 

pjv(ri,r2,...,rAr,r'i,...,r^,f) = \ci{t)\^l[ipi{ri) ^*(r^) + |c2(i)|' [J ^2(ri) (^^(r^) , 

i=l i=l 

where )ci]^ + [02]^ = 1. This function cannot be written in any way as a product of single-particle functions. Hence 

pN is an entangled state, being a kind of the mixed counterpart of the Greenberger-Horne-Zeilinger state [83,84]. 

The principal difference of the two-mode coherent system, wc consider, from the entanglement of trapped atoms, 
studied earlier, is that the nonlinear coherent modes describe not internal states of separate atoms, but collective 
coherent states of the whole system. Therefore this type of entanglement can be called coherent entanglement. 



IX. RELAXATION PROCESS 



In the previous sections, we have considered the situation when Bose-condensed atoms are subject to the permanent 
pumping by an alternating resonant field. One might ask the question, what happens if this resonant pumping is 
stopped? How to describe the behaviour of the system after the coherent resonance? 

If the resonant field ceases acting on the system, say at the moment of time tg, then, after this time, the two- mode 
picture becomes invalid and one needs to return back to the initial equation (6) for the coherent field. This equation 
is to be complemented by the relaxation terms characterizing atomic collisions and possible existence of noncondensed 
atoms [1-3]. Then the coherent-field equation reads 

in^= H[rj] V - y {K2\v\' + KM' + TM) 77 , (106) 

in which H[r]] is the nonlinear Hamiltonian (7), K2 and K3 are two-body and three-body recombination loss-rate 
coefficients [85], F is the loss rate due to the transfer of condensed atoms to noncondensed atoms whose number is M. 
If a cooling mechanism is supported, so that there is an opposite process of transferring atoms from a noncondensed 
cloud to the condensate, then F < is a gain rate. The number of particles in the condensate is iV = \\r]\\-^. The 
initial condition for Eq. (106), in view of relation (8), is prescribed by the function 

r]{r,tQ) = ^ ^{r,to) , (107) 

with (fi{r,t) given by the sum (86) taken at the moment to, when the resonant field stopped acting on the system. 
In addition to loss rates caused by depolarizing collisions, there can exist another internal natural loss rate 

Fi = - I Im (ip, H[p] (108) 

appearing for atoms with attractive interactions, if the number of atoms is larger than the critical number N^. Here 
H[(p] is the nonlinear Hamiltonian (11). Thus, for a spherical trap, the loss rate (108) is estimated [4] as 

Fi ~ 2.867 a;r0(-/V-^c) 



N 
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where N ^ Nc and 6(-) is a unit step function. The loss rates due to two-body and three-body depolarizing collisions 

are 

^ K2 (bp, l^p) , r,^K, |^|3) . (109) 

From the evolution equation (106), we can easily derive the rate equation for condensed atoms, which, being compli- 
mented by the rate equation for noncondensed atoms, makes the general set of equations 

^ = -TiN -r2N'^ -TsN^ -TMN , ^ = -^i{M - Mq) +TNM , (110) 

Civ (liZr 

defining the relaxation process in the system after the coherent resonance. Here 71 is a pumping rate for noncondensed 
atoms and Mq is their stationary number. 

It is not our aim to give in this paper a detailed analysis of the relaxation equations (110). This is a separate problem 
which could be investigated in other publications. Hcrt; our main goal has been to present a thorough description of 
trapped Bose-Einstein condensate under the condition of coherent resonance. Therefore, we limit this section by a 
brief sketch of the way that would allow one to describe in detail what happens after the coherent resonance. The 
relaxation procedure defined by Eqs. (110) can follow quite different patterns depending on the concrete physical 
situation and on the values of the related relaxation rates. In a particular case, when the natural loss rate (108) is zero 
and the are no noncondensed atoms, the relaxation is completely due to depolarizing collisions. Then the relaxation 
time is the lifetime of atoms in a trap. For different traps this time varies between 10 s and 100 s. Note that the loss 
rate (108) is zero for atoms with repulsive interactions and also for atoms with attractive interactions, if the number 
of atoms is less than critical. In such cases, nonlinear coherent modes, created by means of coherent resonance, can 
live, after the pumping resonance field is switched off, quite long time, of the order of the lifetime of atoms in a trap. 
This makes it feasible to study their behaviour as well as to use them for practical applications. 



X. CONCLUSIONS 



We have demonstrated that nonlinear coherent modes of trapped Bose-condensed atoms can be created by means of 
a resonant alternating field. These modes represent collective states corresponding to nonground-state Bose-Einstein 
condensates. Conditions, when such a resonant excitation of coherent modes is possible, are investigated. One 
restriction is that the number of atoms in a mode be less than a limiting number. The latter depends on the type 
of atoms and trap characteristics, and can be as large as 10^. This limiting number is close to that required for the 
stability of atoms with negative scattering lengths. Hence, atoms with attractive interactions can also be employed 
for creating nonlinear coherent modes. Moreover, since the mode limiting number increases for higher modes, it could 
be possible to stabilize a larger number of atoms with negative scattering length by transferring them to such excited 
coherent modes. 

Another restriction on the coherent resonance is that the resonant pumping can last not longer than a critical time, 
before power broadening spoils the resonant picture. However, this limitation is not dangerous because of two reasons. 
First, the critical time for power broadening is rather long, of order 10 — 100 s, which is about the lifetime of atoms 
in a trap. Second, there is no need to pump the system for so long times, since the transfer to an excited coherent 
mode takes essentially shorter times, around 0.1 s. 

Temporal behaviour of fractional populations, in the process of coherent resonance, exhibits dynamic critical phe- 
nomena occurring on a critical line in the parametric manifold. The related time-averaged system displays, on this 
critical line, critical effects typical of phase transitions. The origin of these critical phenomena is the saddle separatrix 
crossing by a starting point of a trajectory. 

Interference patterns and interference current can be observed. These are related to the internal Josephson effect 
and to dynamic barrierless tunneling. 

Atomic squeezing is realized demonstrating that the state coherence is better defined than transition coherence. 
Massive multiparticle entanglement is produced. A specific feature of this coherent entanglement, distinguishing it 
from earlier studies, is that it involves not internal states of individual atoms but collective coherent states of atomic 
condensate. 

After the pumping resonant field is switched off, the behaviour of the system can be described by the rate equations. 
This behaviour depends on the physical setup defining the related relaxation rates. When the loss of atoms in a 
nonlinear coherent mode is caused by depolarizing atomic collisions, the lifetime of the mode is about 10 — 100 s, 
which is the lifetime of atoms in a trap. 

The two-mode picture of a resonant Bose condensate somewhat resembles the two-level picture of a resonant atom. 
However there is a principal diflFerence between these two cases in the nature of the states involved. The notion of a 
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resonant atom involves internal states of a single atom, while that of a resonant condensate has to do with collective 
coherent states of the whole system. Thus, the resonant condensate is a principally new resonant physical system. 
Its similarity with resonant atoms makes it feasible to extend to this type of systems various applications that are 
elaborated for resonant atoms. And its difference from the latter brings hopes of discovering new effects and finding 
novel applications. 
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Appendix 



To illustrate the calculations of Sec. IV, we present here the related explicit expressions for several first nonlinear 
coherent modes, including the ground state mode {n = m = j = 0; p = q = 1), radial dipole mode (n = 1, m = 
j = 0; p = 3, q = 1), basic vortex mode (n = 0, m = 1, j = 0; p = 2, q = 1), and axial dipole mode 
{n = m = 0, j = 1; p = 1, q = 3). The corresponding wave functions are 



V'ioo(r-, </?, z) = ^ "^"°3^°" ^ {uioor^ - l) exp j- i (wioor-^ + vwoz^) | 
'fpoio{r,'P,z) =uoio (-^) ^ ^e"^ exp |- ^ (uqio?'^ + ^^oio-z^) | , 



where w„mj" and Vnmj are the control functions defined by Eq. (39). 
The corresponding integrals Inmj are 



looo = i'^TT)-^^^ = lo = 0M3AM , /ioo = J/o, 410 = ^-^0, Ii 
And the integrals (44) are 



2 -u J -uiu 2 ' ^^^^ ~ ^ ' 



/ 2 , 2 \ / \ 1/2 2 / \ 1/2 

, wioolwo + wfoo) / ■^lOO \ ^ MolO / ■^010 ^ 

^100 = q ; ; , Joio ~ 



JoOl = "001 ^ ^001 ^ 



TT^/^ (Uo +^001) V '^^0 + ^^001 / 

In the weak-coupling limits, when g^/l' — > 0, the control functions behave as 



Mioo ~ 1 - ^ logV^ , uioo - \ logV^ 



Wool - 1 ~ 4 log^ ' ■f^ooi -2^-4 Iog\fi^ ■ 



And for the energy spectrum (42), we have 



£^0 - 1 + 2 + ^o^VJ^ , ^100 - ^ + 2 2 ■ 



iJoio - 2 + ^ + ^ ^offv^ , ^^001 - 1 + Y + ^ iog\/y ■ 
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Then the transition frequency (16), for a transition between the ground state and an excited mode, becomes, respec- 
tively, 

a;ioo,o ^ 2 - 0.031747 g^fv , wqio.o 1 - 0.031747 g^fv , wqoi.o '^v- 0.015874 g^fv . 
And the transition ampHtudes (45) take the values 

ao,ioo — ^0,010 — ^0,001 — , 



aioo.o - aoio.o - 0.031747 g\fv , aooi,o 0.015874 g^fv . 

Hence, in the weak-coupling limit, conditions (46) arc always valid. 
In the strong-coupling limit gv ^ oo, for the control functions we get 

2.282947 2.282947 ^ 

5.823454 1.941151 o 

Uwo — T^T^ , ■^100 — —, 7T7^ ^ ! 



4.565895 2.282947 



2.056114 6.168342 ^ 

which shows that in this Hmit the effective oscillator frequencies diminish. For the energies (43), we find 

2 -I- 18 H- 

£^0 - 0.547539 (gvf^^ + 0.570736 , , E^ioo ^ 0.643949 (gvf'^ + 0.485287 , To/.. . 

(g(i/)2/5 {gvyi^ 



o 1 

Eqw ^ 0.547539 {gvfl^ + 0.570736 ,^ , Eqqi ~ 0.607942 {gvf^ + 0.514029 



2 + 9Z/2 



Notice that passing from the weak-coupling to strong-coupling limit, the effect of level crossing happens, since the 
order of energy levels can be changed. Specifics of the level crossing depend on the value of the aspect ration v. For 
instance, if the trap is cigar-shape {y -C 1) then the arrangements of the energy levels in the weak-coupling limit, 

Eq < Eooi < Eqio < Eioo (5 — * 0) , 

changes, in the strong-coupling limit, to 

Eq < Eoio < Eooi < Eioo {g — > 00) . 

The transition frequencies in the strong-coupling limit are 

wioo.o ^ 0.096410 {gyf/^ , woio.o ^ ' ^001,0 ^ 0.060403 {giyf^^ . 

For the integral (44), we have 

Jioo ^ 0.052069 , Joio = 0.056439 , Jqoi ^ 0.053063 . 
Then the transition amplitudes (45) become 

ao,ioo ^ 0.140202 {gyf/^ , aioo,o ^ 0.101636 {guf/^ , 
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ao,oio ^ 0.170345 {guf^ , aoio.o ^ 0.170345 {guf/^ , 



ao,ooi - 0.147059 {giyf/^ , aooi,o - 0.122898 (giy)^/^ . 

This shows that in the asymptotic hmit gv oo, inequahties (46) do not hold. They become vahd outside the region 
of convergence of this asymptotic Hmit. The corresponding condition (48), for the considered lower modes, yields 

\g^\<l.l + {p=l,q = l), 

\gu\<0.7 {iS + u^f^ (p = 3,q=l), 
|H<1-1 + (p = 2, g = l), 

|.gz/| < 0.8 (2 + 9i/2)^/^ (p=l^q = ^). 

These inequalities demonstrate that large values of the coupling parameter g I can always be compensated by 
accepting an appropriate trap shape, being either cigar-shaped {u <^ 1) or pancake-shaped {u 1). 
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Figure Captions 



Fig. 1. Characteristic evolution of the phase portrait on the population difference, s{t), - phase difference, x{t), 
plane, under zero detuning £ = and varying amplitude of the resonant field: (a) b = 0.1, (b) b = 0.4; (c) b = 0.49; 
(d) b = 0.51; (e) b = 0.8. 



Fig. 2. Dramatic change in the temporal behaviour of the population difference s{t) (dashed line) and phase 
difference x{t) (solid line) when crossing the critical line at the point {be = 0.5, £ = 0}. The initial conditions are 
So = —1, xo = 0. Time is presented in dimensionless units, as explained in the text. In Figs. 2(a) and 2(b), the value 
of b is, respectively, just below and just above the critical point be, with the variation of Ab = 10~^. 



Fig. 3. The change of dynamics in the population difference s{t) (dashed line) and phase difference x{t) (solid 
line) when crossing the critical line at the point {be = —0.5, e = 0}. The initial conditions are sq = —1, xq = 0. Time 
is dimensionless (see the text). In Figs. 3(a) and 3(b), the value of b is, respectively, slightly above and slightly below 
the critical point be, with the variation A6 ~ 10~^. 



Fig. 4. The squeezing factor as a function of time for the transition amplitude b = 0.49, zero detuning, and 
the initial conditions sq = —1 and a;o = 0. 
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